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Abstract 

A method is presented for analysis of correlation function of two non-identical particles with 
strong and Coulomb interactions, resonance formation, channel coupling and spin structure. For 
resonance reactions we derive a formula giving the small distance contribution to the correlation 
function. The formalism is used to analyze the preliminary RHIC data on ir^r^ correlation 
measurements. The H*(1530) resonance is successfully described. The 7rH source size is obtained. 

PACS numbers: 25.75.-q, 25.75.Gz 
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I. INTRODUCTION 



Measurements of momentum correlations of two low relative momentum particles pro- 
duced in heavy ion collisions provide a unique information on spatio-temporal picture of the 
emission source at the level of few fermis. The study of the collision region on the femtome- 
ter scale via two particles correlations is called femtoscopy (see, for example, [^-[3]). At 
the early stages the studies were focused on the identical pions correlations arising from the 
wave function symmetric g Th.s type of analysis has a deep ana.ogy with the «HBT 
interferometry" used in astronomy [9j. Since then, measurements have been performed for 
different systems of both identical and non-identical hadrons. High-statistics data sets were 
accumulated in heavy ion experiments at AGS, SPS and RHIC accelerators [lO]- Corre- 
lations are significantly affected by the Coulomb and/or strong final state interaction (FSI) 
between outgoing particles. The non-identical particle correlations due to the FSI provide 
information not only about space-time characteristics of the emitting source, but also about 
the average relative space-time separation between the emission points of the two particle 
species in the pair rest frame (PRF) 18]. 



Maybe the most exotic system studied recently by STAR collaboration is 7r — S 19, 



20j: the particles composing the pair have one order of magnitude difference in mass plus 
AB=1/AS=2 gap in baryon/strangeness quantum numbers. It is challenging to study FSI 
of such exotic meson-baryon system and to extract information about the 7r — H scattering 
lengths. The other important reason to study 7r — S correlations is that multistrange baryons 
are expected to decouple earlier, than other particle species because of their small hadronic 



cross-sections 



2l| . allowing one to extract the space-time interval between the different 



stages of the fireball evolution. 



Preliminary results for the 7rS system are available from STAR Collaboration 19|, |20 ] . 
The following important observations were made: 

• Decomposition of the correlation function C(k) = C(k, cos 8, ip) (k is the center-of- 
mass (cms) momentum), from 10% of the most central Au+Au collisions into spherical 
harmonics, provided the first preliminary values of R = (6.7 ± 1.0) fm and A out = 
(—5.6 ± 1.0) fm. The negative value of the shift parameter A out indicates that the 
average emission point of H is positioned more to the outside of the fireball than the 
average emission point of the pion. 
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In addition to Coulomb interaction present in previous non-identical particle analysis 
the 7r + S~ correlations at small relative momenta provide sufficiently clear signal of the 
strong FSI that reveals itself in a peak corresponding to the H*(1530) resonance. The 
peak's centrality dependence shows a high sensitivity to the source size. 
In 



19 



20j was shown qualitative agreement with the model calculations in the Coulomb 



region and overestimation of the peak in the S*(1530)-region. 

It is clear that the problem of 7rS correlations deserves a thorough theoretical treatment 
which is undertaken in the present paper. We will be concentrated on finding the expression 
for the 7rS wave function (WF) with Coulomb and strong interactions included. 

The wave function enters as a building block into a correlation function (CF): 

C(k) = J dfSir)^^ 2 . (1) 

here S(r) is a source function, ^jg(f) is the two final-state particles WF. 

The influence of the Coulomb interaction can be taken into account following the standard 
procedure described in textbooks. A model independent approach to strong FSI is absent 
since the low energy hadron interactions can not be described from the first principles of 
QCD. Phenomenological approach to the combined treatment of Coulomb and strong FSI 
is based on the effective range expansion of the strong amplitude. In case of 7rS system the 



problem is rather complicated due to the following factors: 22l|-(25|. 

• The superposition of strong and Coulomb interactions 

• The presence of H*(1530) resonance 

• The spin structure of the w.f. including spin-flip. 

• The fact that the 7r + H~ state is a superposition of I = 1/2 and 1 = 3/2 isospin states 
and that 7r + 2~ state is coupled to the 7r°S° and that the thresholds of the two channels 
are non-degenerate. 

• The contribution from inner potential region where the structure of the strong inter- 
action is unknown. 

The description of the 7rS correlation is a twofold problem. First of all one has to 
construct the nE WF with all factors enumerated above included. Secondly, one has to use 
a reasonable model for the source function. In the present paper we concentrate on the 
first task, while the source function is taken in the simplest Gaussian form. Limitations of 



the naive Gaussian model are well known and will be briefly discussed at the end of the 



paper. In 
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FSI model from 



the Blast Wave model for the source was used in combination with the 



221 ]. the simultaneous description of the Coulomb and resonance regions 



was not obtained. Here we obtained successful description of the source size using simple 
Gaussian source model. It is clear that more elaborated source model is need to describe 
simultaneously source sizes and emission asymmetries. 

The paper is organized as follows. Section |TT] is devoted to a somewhat pedagogical 
introduction to the description of the FSI. Our purpose will be to expose the formalism 
in a way suited for the construction of the CF. In Section II III we consider the spin-isospin 
structure of the WF and the coupling between charged and neutral 7rS channels. Here we also 
present the first fit to the experimental data. This fit encounters problem in the resonance 
region, and the blame for this discrepancy lies in the disregard of the small distance region. 
The small distance contribution to the WF is derived in Section [TV] Section [V] is devoted 
to the analysis of the experimental data. In Section [VTl the main results are presented and 
open problems are formulated. 



II. THE BASIC FSI FORMALISM. 



As it was adumbrated above, the structure of the 7rS WF is rather complicated. To set 
the scene for the detailed treatment of the FSI in the 7rS system we start with a pedagogical 
introduction into the FSI problem in the context of femtoscopy. Almost all results presented 
in this paragraph may be found in the literature. What we have tried here is to gather in 
one place the properties of the out-state WF \f) needed to calculate the CF. Thus our 
treatment does not seriously overlap with any other in the literature. 

There are two complete sets ty^\r) of the continuous spectrum wave functions. In 
-(f) the ingoing flux is nonzero only in the direction n = —rh, n = k/k, m = f/r. This 
corresponds to particles moving along n toward the center. The flux of outgoing particles in 
vj/v") (f) consists of two parts: non-scattered ones moving along n in the direction n = m out 
of the center, and scattered particles moving in all directions. On the contrary, in (r) 
ingoing particles are moving from all directions and outgoing particles are propagating only 
along the ray n = in. Asymptotically at r — > oo both functions (r) contain the same 
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plane wave e tkr : 

ikr 

(f) ~ e iW + f{nm) — (2) 

*i _) (r) ~ f(-nm)— (3) 

Coulomb interaction and channel coupling bring some distortions into these equations - see 
below. The following relation between ^^\r) and ^ ^(r) holds 

rH/,* _ ,T,(+) 



*r(*o = *T(f) (4) 



The wave function (f) which is called the out-state is used to describe particles 
produced in some process. Suppose, for example, that some state <f>(f) is created in the 
potential U(f). Then the flow of the particles with the wave vector k emitted from the 
center is given by 

N(k) = J dr(j>{r)^\r)\*. (5) 

The out-state \r) has to be used in Eq. ([1]) for the correlation function. 

Consider a pair of particles interacting via attractive Coulomb potential (like tt + — S~). 
The WF ^"Vf) reads 

kc 

Vt\r) = e^ /2 r(l + iriy^Fi-ir], 1, -i{kr + kr)). (6) 

Here rj = l/kas, clb is the Bohr radius, a# = 1/a/x (=214 fm for 7r + — fi is the 
reduced mass (=126 MeV for n + — Br), a = 1/137, F(a,j,z) is the confluent hyper- 
geometric function. The region rj > 1 is called the atomic energy range. In this region 
Coulomb interaction dominates even in presence of strong interaction. The wave function 
is normalized so that 

where c 2 (k) is Gamow's factor. 

As an illustration let us take in ([T]) a Gaussian model for the source function 

S(r) = (STr^R'y 1 exp(-r 2 /4R 2 ). (8) 

Expanding ^j- (f) in r/as and keeping only the first order term we obtain the following 
result for the correlation function (pQ): 

C{k) = c\k){l-^ + ^}. (9) 
V7T a B 



The partial wave expansion of the Coulomb WF is given by 

= VrP Q il{2l + Ve-^lM* ( Yr) " ^ 
where ai is the Coulomb phase 



ai = argY{l + 1 — if)), 



-to 'i 



F(l+l+iri) 



\r(i+i+i V )\ ' 

and Fi(r), kr) is the regular solution of the Coulomb problem 



Fjfo kr) = N ie - ikr (kr) l+1 F(l + l + ir],2l + 2, 2ikr), (12) 



,i 1 



(11) 



Ni ~{^)II {m' + v 2 )} 112 . (13) 

^ ' m=l 

Asymptotically F(rj, kr) behaves as 

Fi(i], kr) ~ sin ^kr - y + r]ln(2kr) + 07^ . (14) 
Irregular solution is denoted by Gi(rj, kr) and it has the following asymptotes 

Gi(rj, kr) ~ cos ^kr — — + r)ln(2kr) + a^j . (15) 

Asymptotically at £ = r(l + cos#) 3> 1/fc one has 

(?? 2 \ - fc~\6) 
1 + i— J expjifcr + ir]ln(k^)} H — exp{— iAr — ir]ln(2kr)} , (16) 



^'w = 2tcL»|rS-S ° p{ - iii, ^ (c "l )} - (17) 

When the produced particles interact only via strong forces their wave function outside 
the interaction range (at r > e ~ 1 fm ) has the following form 

^\r) = yX/ {21 + • (18) 

The radial wave function Xki(j') asymptotically behaves as 

Xki(r) ~ sin (kr - y + 5 Z j (19) 



Consider now a situation when FSI is caused by the combined action of the Coulomb and 
strong interactions. Then at r > e we have 



= Vr E^ 2/ + f T~ ) • (20) 



kr ^— ' \ kr 

1=0 



The radial WF Ukiir) is a regular solution of the Schrodinger equation containing the sum 
of Coulomb and strong potentials. An important point is that 5i entering into (120]) is not 
identical with the pure hadronic phase shift 5i in (fl8|) . The difference is called the Coulomb 



correction 



26] and will be explicitly introduced later. 



At r > e the radial wave function in ( 1201) has the form 

u k i{r) = Fi(r), kr) cos5 z + Gi(r), kr) sin 5/. (21) 
Making use of f l2Tl) and of the expression ( flOl) for the Coulomb WF we obtain 



= *&(^ + 3SF E£o<'(H + l)e~^(e-^ - l)[F,fo, fcr) + <Gf,(i|, A;r)]P, (f£ 

= + ESo JgviikrmW (g) , (22) 



where 



, , . v^47r , ■ f Fi + iGi \ 
<Pi(kr) = ^—% l e- lai , , (23) 



kr \ % 

p -2iSi _ 1 

Ti(k) = (24) 



Recalling the asymptotic form of ^(r) given by Eqs. (|T6|) - (|T7|) and the asymptotic behavior 



fee 



of -F) and defined by Eqs. Q14p -( jl5p one arrives at the following asymptotic of ^1 (f 



fees 



~ ( X + ^) e ifc>+ ^ n(fcf) + (ft\0) + f£\9)) ^ e - i{kr+r,ln2kr \ (25) 

where the Coulomb amplitude fc (9) is given by (IT7I) and the Coulomb modified strong 
amplitude fes (9) is equal to 

] (9) = ^ f> 1} ' (2/ + 1 ) c ~ ato, ( c ~^ - W f^J • ( 26 ) 



Now let us turn to the low energy expansion of the WF (1221) . For the moment we are not 
interested in its spin-isospin structure. This problem will be addressed in the next Section. 
Modification of the effective range expansion due to Coulomb interaction is a problem whose 



29| . Here we quote the results relevant for the construction of the 



solution is well known [27 
out-state WF ^ \r). Coulomb corrections are most important for the S- wave. The basic 
role in the low energy expansion is played by the effective range function M(k). In S"-wave 
we have 

kc 2 {k) 



-2iS 



To 



1 



-2« Mfjk) + A + ik 
Here A is the so-called Schwinger's correction |30[ 

A = — In — + 27 



a B 



a B 



(27) 



(28) 



7 ~ 0.577 is Eugler's constant, k is the Coulomb corrected momentum. 

2i 



k = c\k)k 



a B 



(29) 



(30) 



h(x) = ln(x) + Re^ ( 1 + % 



•'' 1 I -7 + ln(x) + ^ x»l 

with ^f(z) = 4-lnT(z). In what follows the Schwinger's correction will not be taken into 
account. The function M(k) has the usual effective range expansion 

M(k) = ^ + ^R Q k 2 + .... (31) 

Similar procedure for the P-wave leads to the result 



e 



-2iS 



1 jfe 3 (l +7] 2 )d 2 (k) 



with 



T = r : = ^ J >- ^ (32) 

-2i AT (A;) + ikk 2 (l + r^ 2 ) 



iV(A;) = ^ + ^i?iA; 2 + . . . . (33) 



where b is the scattering volume with the dimension fm 3 , and Ri = —3/d with d having 
the interpretation of the range of forces. The Coulomb correction in P-wave similar to A 
(see ( 126]) ) can be neglected. 

The important feature of the 7rS system is the existence of the P-wave resonance 
S(1530)Pi3 with a width T = 9.1 MeV. The corresponding value of the cm. momentum 
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is k = 150 MeV, rj = 1/kas — 1/160 and therefore Coulomb corrections in this region are 
small. Let us recast the P-wave term fl32|) into the Breit-Wigner form. We get 

r k\l + v -)c\k) ^ 1 r 

N(k) +ikk 2 (l + V 2 ) 2E-ER-IT/2' 

where 

e = A(s -^-'" 2 > „ 2„E, (35) 
E = </s-(M + m), (36) 
E R = M* - (M + m), (37) 

r = ,38) 

with M and m being the masses of H and 7r correspondingly, M* = 1530 MeV is the mass 
of the S* resonance, /i = mM(M + m) . The scattering volume b is expressed in terms of 
the Breit-Wigner resonance as 

" ^ ^ (39) 

III. SPIN-ISOSPIN STRUCTURE OF THE WAVE FUNCTION. 

Up to now the spin-isospin structure of the WF has been disregarded. In this 

section we shall wind up this lacuna. 

It is sensible to begin by considering the angular momentum decomposition of the WF 
with a given isospin /. The isospin structure of the WF as well as the coupling to the 7r°S° 
channel will be included at the next stage. The amplitude of meson-baryon scattering has 
the following form: 

T = f + g(an) (40) 

n = (41) 
I [Mil 

The second term in ( I40p corresponds to the spin-flip amplitude 1 . Now let us write down 
the expression for the WF with the spin-flip term included. This expression will replace 



[1] The importance of the elastic spin-flip and charge-exchange amplitudes in this problem was indicated to 
the authors by R. Lednicky. 
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the spin-less equations (!20|) . fl22|) . For any type of the interaction the conserved quantum 
numbers are j 2 , I 2 and j z = rrij. The WF having definite values of the above quantum 
numbers is constructed from the direct product of the WFs (I2"0l) . (l2"2"j) and the spin WFs 











H 

















(42) 



corresponding to rrij = m s = +1/2 and rrij = m s = —1/2 respectively. The resulting 
expressions read 



= *<.-' I 1 1 + 







— / VT+1Y 1>0 \ _ r-r / VlY l>0 



+ I I ' — I + I I I (43) 




1 



+ e,^m I _ r rr , , . (44) 



Here the symbols (Z±) correspond to j = I ± 1/2. In deriving (14111) -( 1441) use was made of 
the relation Y* m (-k/k) = (-l) l+m Y lt _ m (k/k). 

The low energy region of 7rS interaction up to the S*(1530) resonance is dominated by 
S- and P- waves. Keeping only these two amplitudes we can rewrite (1431) as follows 




(45) 



Similar expression holds for rrij = —1/2. The last term in fT45]) corresponds to spin-flip 
and this contribution vanishes if 5i+ = 5i-. Recalling about isospin doubling we conclude 
that the WF fj4"5"]) contains two S-wave amplitudes and four P-wave ones. To determine 
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six amplitudes from femtoscopic experiments is hardly possible. To reduce the number of 
parameters we shall assume that the dominant interaction in P-wave occurs in a state with 
j = I + 1/2 = 3/2,1 = 1/2 containing the S*(1530) resonance (the amplitude T 1+ in (H5T) ). 
Then the number of parameters becomes equal to three. The actual number is two since the 
parameters of the S*(1530) resonance are known from the experiment. There might be one 
more parameter. The point is that expressions (1201) . (1221 . (I43H45I) for the WF correspond 
to distances r > e ~ 1 fm where the strong interaction is assumed to vanish. We shall 
address this problem in the next section and show that the contribution of the inner region 
is important while e is not a relevant parameter. 

We now proceed to discuss the isospin structure of the WF and coupling between charged 
and neutral 7rS channels. There are two 7rS systems with opposite electric charges: 7r + 2~ 
and 7T~2~. For defmitiveness we shall consider the charged channel 7r + S~ which is coupled to 
the neutral channel 7r°S°. The corresponding thresholds are not degenerate with the 7r + 2~ 
threshold being by ~ 11 MeV higher. The problem can be treated either in the channel 
basis |7r°S )}, or in the isospin basis {\I = |), \I = §)}. The relation between the 

two frames is given by 





l ) \ * / 11 = l ) 

2/ I _ TT I 1 2/ 



, U\ ' 11 . (46) 

3 W = 



by several authors 28] 



Shaw and Ross 



28 



Strong interaction is diagonal in the isospin basis while we are interested in the channel 
WF |/T + 2~) (we shall use the subscript a for the quantities corresponding to this channel 
and P — for shose corresponding to |7r°H )). Similar problems have been treated previously 

29], 3JJ. We shall more or less follow the approach proposed by 
]. In Section [TT] we have introduced the function M(k) = k cot 5 which 
allows to perform the effective range expansion in presence of the Coulomb interaction (see 
Eqs.( l2~TH3~T|) ). The generalization to the two-channel system reads 

f = k l+1 /\M - ik 2l+l Y x k l ^\ (47) 

where 

JU(*"°V M=( M " M A (48) 
\0 kg J \ M ea Mpp J 

Consider first the S'-wave. In single channel case in the scattering length approximation one 



11 



has M = 1/a (see |3"I1) . In the two-channel case in the isospin basis one has 

M 1 = \ ~ s ° \ , (49) 




where a s corresponds to I = 1/2, and at — to I = 3/2. Transformation to the channel basis 
reads 

M = UM I U- 1 = -\ _; s a * V a ° a *J I (50) 

3 \ y/2 -- + - - + - / 

where the matrix U was introduced in (|46j) . According to (|47[) the two-channel generalization 
of (27) is 

T = k a c 2 (k a )(M^ + tk l3 )d-\ (51) 
d = (M aa + ik a ){Mpp + ikp) - \i nS .\l -„-,. (52) 
where as before the Schwinger's correction ( 128]) has been neglected. We can now substitute 



(joTj) into fj45|) and obtain the closed expression for the S'-wave component of the WF. How- 
ever, in the two-channel case this would not be a complete answer. One should add the 



contribution from the coupling to the |7r°S°) cnannel, i.e., a term *H ~ SV a ¥ 8 I The 



corresponding contribution to (145]) is 

1/2 



<Li/ 2 (^ <" «)s = ^oo^o r • (53) 



Here 



^(V/?) = ^f' e -^(^)/ i H(^) ) (54) 

KpT p 

Ro = -c(k a )(k a kp)V 2 M a pd-\ (55) 
where hi(z) is Hankel function 

hf\z) = e ±iz , h^\z) = (l±i/z)e ±i{z - n/2) (56) 



We note that ooi(p) = lim^o fi{p)i where ifi(p) was introduced by 

As it was stated above the .S'-wave component of the |7r + S~) WF depends upon two 
parameters, namely the isospin scattering lengths a s and a t . With our sign convention 
moderate attraction corresponds to positive signs of the scattering lengths. Next we turn 
to the P-wave. We assume that the S*(1530) resonance with 7=1/2 plays the dominant 
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role and therefore we neglect the P-wave with I = 3/2. The resonance is coupled to both 
7r + 2~ and 7r°H° channels. The corresponding amplitude T\ + and R\ + are evaluated using 
the matrix U (1461) . Denoting the resonance width in the isospin basis by V we make the 
transition to the channel basis 



(57) 



In line with (|34|) we shall omit the Coulomb corrections and write 

= _i (uru-% = _i r (58) 

2E-E R -i\ 3E-E R -q< 




ft + = "fr = ^ „ ,r - (59) 



2E-E R -%\ ^E-E R -i\ 
In writing ( |59|) we made the approximation 



^U^r^y 2 ^) 1 '^-^. (60, 

Two remarks are in order at this point. First, the amplitudes Ti and Ri are added non- 
coherently. This will be visualized in the expression for l^^^ 2 . Second, in order to 
calculate the CF according with |T|, or a similar equation, one has to take rp = r a = r. 
This is a reasonable approximation since the reduced masses in 7r + H~ and 7r°S° channels 
are close to each other. 

Now we can collect all pieces together and write the expressions for ^m] and | | 2 with 
Coulomb interaction included in all partial waves and strong interaction in S and P-waves. 
It is easy to see that 

Itf(-) |2 = hJrH |2 ( 61 ) 
I mj=l/2l I irij=— 1/21 V / 



and therefore 



So, we write 



* h i 2 =~ £ i*yi 2 =i<u /2 i 2 ( 62 ) 

m 3 =±l/2 




<U = M n )+^| n )+VsM v, . ) + VsM Yl _, I- (63) 
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|^(-)| 2 = Y 00 (K| 2 + |^| 2 ) + -^Y 00 Y 10 {ReA a B* a + ReApBp) + 

+*y* (\B a \* + W) + f l^wi 2 (i^i 2 + i^i 2 ) = 

= i (\A a \ 2 + \A?\ 2 + |P Q | 2 + l^l 2 ) + (ReA a B* a + ReApB*) + 

+^(|P Q | 2 +|P /3 | 2 ). (64) 

Here 

A a = + VoT , B a = Vl T 1+ , (65) 

/ , \ 1/2 / , x 1/2 

A, = ( c^Po, 5/3 = f f\ wii2i+. (66) 

Equation ( 1S"4"1) is a sought for result which should allow to evaluate the CF according to (00). 
In the next paragraph we shall see that this is still not a complete answer. 



IV. INCORPORATING THE INNER REGION. 



With the WF f l64j) at hand we can use Eq.([T]) and try to fit the experimental data 
on the CF. The results are shown in Fig. [I] by dashed line for the fireball radii 
R = 7.0 fm and zero scattering lengths. The low momentum region is properly described. 
However, in the region of the 5* (1530) resonance we see a dip-bump structure instead of 
the experimentally observed resonant behavior. The wiggly behavior of the calculated CF 
is explained by the interference at 6 ~ n between ¥^ and the resonant P-wave, see (1221) . 
or (T441) 2 . Coulomb effects are not important here and for the sake of clarity they can be 
neglected, the phenomenon is present already in the general expression ([3]). Keeping in ( 122]) 
the resonant P-wave one has 

V>(->(0 = 7r) ~e" ifcr (1 + t 1 - - -]. (67) 

y v ; V krE R -E + iT/2J K ' 

This leads to the interference term 

i / (-) 12 2 r(p i? - e) 

^ ~ MP-p R ) 2 + r 2 /4 - (68) 



[2] The interference at ~ 7r between the incident wave and the scattred one was observed and explained in 
22|. 
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The last expression clearly displays the dip-bump behavior seen in our calculations. The 
same result can be derived from the expression PI) for |^ ( - } | 2 . 

Contrary to this result the experimental points shown in Fig. [TJ exhibit a resonant struc- 
ture. It means that something is missing in our calculations. Recall that the WF-s ( 14"3"|) , 
OS]), (JM]) correspond to the region r > e where the strong potential is assumed to vanish. The 
contribution from the region r < e is proportional to the time the impinging wave packet 
spends there. For the resonance this time is much larger than is needed to cross the sphere 
with the radius r = e. 

We come now to a formal account of the small distance contribution. First, we split the 
CF ([T]) into two parts 

C(k)~ I dfS(r)\^t\r)\ 2 + S(0) j dr\¥^\r}\ 2 . (69) 

J r>t J r<e 

The source function was taken out of the second integral because S(r) has a characteristic 
scale R > 4fm ^> e. The integral over the region r < e entering into the last equation can 
not be evaluated directly since the structure of the WF at small distances is unknown. 



However this integral can 
Luders-Wigner formula 



:>e expressed in terms of the scattering phase shifts using the 
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3 lM3 31] . This approach is simple and effective when FS 



is 



dominated by a narrow resonance. Otherwise one can apply the method proposed in [24] 
(see also |22j). 

Consider the state with j = 3/2, 1=1/2 containing the 2*(1530) resonance. For the 
radial WF Xj,i( r ) — X3/2,i/2( r ) the following relation holds 

J = [ ^IX3/ 2 ,i/ 2 (r)| 2 = ± + e - i- sin2(fce + 6 3/2 , 1/2 - tt/2)) . (70) 



The phase is given by 



53/2,1/2 = ~ arct9 2(E - E R ) ' ^ 



From (1?0D and (1711) we obtain after a short calculation 



o. 2 j _ k r / 2 | ({E - E R f - T 2 /4) sin 2ke - T(E - E R ) cos ke 

fi(E — R R ) 2 + T 2 /4 + 2k ( (E - E R y + rV4) 1 ) 

where /x is the 7rS reduced mass. This expression may be to a high accuracy approximated 
by the first term only. Indeed, at E = E R the sum of the last two terms yields 

sin2ke 2 2 



2k 3 
15 



(ke) 2 e ~ e, (73) 



since (ke) ~ 1. The first term at E = Er is given by 

2k 2 

^- f » e . (74) 

Physically (I74p means that particles spend inside the small distance region the time ~ 1/T 
which is much larger than the time e/v = efi/k needed to cross this region. 

Therefore from (I70p and ( I64p we obtain for the inner region contribution entering into 
( |69l) the following result 

UdA^-\r)\* = \U d ^W^\ 2 = U dQ ( 1± ^) Jo dr(4n\ Xw (r) | 2 ) = 

= * / e ^IX3/2,i/2(r)| 2 = fe pz^Ti - (75) 

Here the factor 2/3 is the isospin Clebsch-Gordan coefficient squared, the factor Air in front 
°f |X3/2, i/2( r ) | 2 results from our normalization condition, see (I2"3"j) and (T54"j) . We can rewrite 
this relation in a relativistic form, as 

df|tf(-)(f1l 2 = — — C76) 



s 2 - (M 2 -mW 



^rel „2 ™2 , \2> (77) 

with M and m being the S~ and 7r + masses respectively, out of the resonance region (s+M R ) 
in (ITT)) is replaced by 2s. 

In Fig. [T] by solid line we show the CF calculated according to ( |69l) with the inner region 
included. 

V. ANALYSIS OF THE EXPERIMENTAL DATA. 

We now turn to the comparison of the calculated CF with the experimental data presented 

in 
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20J. From Fig. [TJwe conclude that the data are fairly well described if we take the 



fireball radius R = 7.7 ± 0.7 fm. The fitting was performed by the minimization of the 
functional \ 2 — ^2(CF(k) exp — C F(k) mo d e i) 2 / &l . The only fitted parameter was R, the 
scattering lengths a s and a t were put equal to zero (see below), calculations were performed 
for e = 0.8 fm, but as shown above for a narrow resonance the results are stable with respect 
to variations of e at e <C R- 

Interpreting Fig. [TJ one should keep in mind that in different momentum intervals the CF 
yields different physical information. Unless 5-wave scattering lengths are large, which can 
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hardly be expected for the 7rS system, and unless the fireball radius is anomalously small, 
attractive Coulomb interaction becomes dominant already at k < 0.05 GeV/c. Strong in- 
teraction effects in this region are roughly speaking proportional to af/R 2 (i = s,t) and 
therefore negligible for Oj <C R. Therefore a good fit in Fig. [T] with R = 7.7 ± 0.7 fm is 
obtained with zero scattering lengths. Only very precise experiments may allow to deter- 
mine the values of a s and at- This conclusion is illustrated by Fig. [2] which shows the low 
momentum CF for different values of R and a s , at 3 . Given the present accuracy of the data 
we simply put the scattering lengths equal to zero and ignored correlations of their errors 
with errors of the source size. 

Important information on the fireball radius and on the FSI can be inferred from the 
resonance region. An important point is that the resonance phenomena is very sensitive 
to the source size. This is illustrated by Fig. [3j The fit to the resonance effect allows to 
obtain the source size value rather reliably. Based on the data points presented in It], [^(J 
we conclude that the source size is R = 7.7 ± 0.7 fm which agrees with the value obtained 



m 



19 



2CH ] where the low-k , Coulomb dominated part of the CF was selected for fitting, 



excluding the region of the peak. Unlike |19|, |20|] we succeeded in describing H*(1530) 
resonance region. 

It should be however noted that the source model which we used here (Gaussian in PRF) 
is oversimplified. The first step toward a more realible model of the source function is 
suggested by the expression (!6i!) for the WF squared. This expression is tantamount to 
the expansion in terms of spherical harmonics with I = 0,1,2. Therefore the CF can be 
decomposed with the I = 1 term giving the emission point difference between 7r and H. Such 
analysis was done in [3]. Other important points left beyond our present investigation are 
the influence of pions produced from resonance decays and non-Gaussian tail of the source 
unction. Both factors are significant for pion and kaon production in heavy ion collisions 
35)]- {37(1. Experimental data on 7rS correlations were analyzed in {19II20I1 with the account of 
these effects. Roughly speaking they lead to the increase of the source size but the problem 
has to be carefully investigated. The same is true for the distortion of the CF by the collective 



flows. As demonstrated in 



for d — a correlations with resonance formation collective 



flow may reduce the source size and lower the resonance peak. One should, however, keep 



[3] Recent lattice calculations yield the result a(7r+2°) = -0.098 ± 0.017 fm [34| 
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in mind that direct production of S*(1530) may compensate this effect. 
VI. SUMMARY. 

We have developed a formalism to evaluate the CF for the system which comprises several 
features aggravating the analysis. These features are: superposition of strong and Coulomb 
interaction, presence of a resonance, spin structure, channel coupling. We have also pro- 
posed a new method to incorporate the contribution from the small distances where the 
structure of the interaction is unknown. This technique which is legitimate for large enough 
source size has been successfully applied to the preliminary RHIC data on 7r + 2~ correla- 
tion measurements. The Gaussian source in the pair reference frame was considered. The 
resulting value of the source size radius is equal to R = 7.7 ±0.7 fm. There is no sensitivity 
of the data to the values of the 7r + 2~ S-wave scattering lengths. There is a large sensitivity 
of the correlation function to the source size in the region of S*(1530) resonance. A good 
description of the S*(1530) resonance region was obtained. The developed formalism can 
be useful for femtoscopy studies of the systems with the narrow resonances born due to FSI 
e.g. K+K- (K + K~ ->■ ). 
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FIG. 1: The CF of 7r + H~ system as a function of PRF momentum k for R = 7.7 fm and 
zero scattering lengths without including the inner region correction (dashed line) and the full 
calculations (solid line), the STAR collaboration experimental data points from 19|, 
stars). 
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